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>> ' Abstract 

We consider lattice field theories with topological actions, which are 
invariant against small deformations of the fields. Some of these actions 
have infinite barriers separating different topological sectors. Topolog- 
q ' ical actions do not have the correct classical continuum limit and they 

cannot be treated using perturbation theory, but they still yield the cor- 
rect quantum continuum limit. To show this, we present analytic studies 
of the 1-d 0(2) and 0(3) model, as well as Monte Carlo simulations of 
the 2-d 0(3) model using topological lattice actions. Some topological 
actions obey and others violate a lattice Schwarz inequality between the 
action and the topological charge Q. Irrespective of this, in the 2-d 
0(3) model the topological susceptibility Xt — (Q 2 ) /V is logarithmi- 
cally divergent in the continuum limit. Still, at non-zero distance the 
correlator of the topological charge density has a finite continuum limit 
which is consistent with analytic predictions. Our study shows explicitly 
that some classically important features of an action are irrelevant for 
reaching the correct quantum continuum limit. 



"Dedicated to Ferenc Niedermayer on the occasion of his 65th birthday 
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1 Introduction 



Universality is a key concept in classical statistical mechanics and in quantum field 
theory. In particular, in lattice field theory numerous lattice actions yield the same 
universal continuum limit. It is well known that locality is vital for the viability 
of universality. A universality class is characterized by the space-time dimension 
and the symmetries of the relevant fields. In order to construct lattice theories that 
fall into a desired universality class, one often imposes additional features of the 
corresponding classical theory on the lattice action. Usually one constructs a lattice 
action by replacing derivatives of the continuum fields by finite differences of the 
lattice fields. Such a discretization ensures the correct classical continuum limit. 
In addition, the lattice theory can then be investigated using perturbation theory. 
For example, for QCD it has been proved rigorously that the lattice regularization 
yields the same continuum limit as perturbative regularization schemes, such as 
dimensional regularization [l|, 0]. Taking advantage of universality, and following 
Symanzik's improvement program 0,0], one can systematically construct improved 
lattice actions p, 0| which eliminate lattice artifacts up to a given order in the lattice 
spacing. At a fixed point of the renormalization group, so-called classically perfect 
lattice actions have been constructed, which are free of lattice artifacts at the clas- 
sical level p7l4lo| . In particular, in asymptotically free theories, including QCD and 
the 2-d 0(3) model, a classically perfect fixed point action has been constructed 
by solving a minimization problem. In this paper, we proceed in a very different 
direction. In order to test the robustness of universality, we explicitly violate classi- 
cally important properties of the action, such as the classical continuum limit, the 
applicability of perturbation theory, or the Schwarz inequality between the action 
and the topological charge. Hence, in contrast to Symanizik's lattice actions, which 
may, for example, be 1-loop improved, the actions that we will study can be viewed 
as tree-level impaired, but they are certainly still local. As we will see, even with- 
out appropriate classical features, the lattice theory acquires the correct quantum 
continuum limit. This also holds in one dimension (i.e. in quantum mechanics), al- 
though in this case one does usually not rely on universality. From the point of view 
of Wilson's renormalization group applied to critical phenomena the irrelevance of 
classical features of a given lattice action is perhaps not too surprising. For example, 
it is well known that the 4-d Ising model, whose classical Hamiltonian does not have 
a meaningful continuum limit, is in the same universality class as the 4> 4 quantum 
field theory. 

We will investigate local lattice actions that are invariant against small continu- 
ous deformations of the lattice fields. Such actions — which we will call topological 
lattice actions — have infinitely many flat directions because they do not suppress 
small field fluctuations. As a consequence, they do not have the correct classical 
continuum limit and perturbation theory is not applicable. Depending on the nature 
of a topological action, it may or may not obey a Schwarz inequality. In this paper, 
we study two different types of topological lattice actions. The first one constrains 
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the angle <p between nearest-neighbor O(N) spins to \(p\ < 5. All field configurations 
that satisfy this constraint have the same action value 5 = 0. Besides not having the 
correct classical continuum limit, this lattice action violates the Schwarz inequality 
between action and topological charge. The quantum continuum limit is reached by 
sending the maximally allowed angle 8 to zero. Patrascioiu and Seiler as well as 
Aizenman li[ have used an action with an angle constraint to simplify the proof of 



the existence of a massless phase in the 2-d 0(2) model. Furthermore, Patrascioiu 
and Seiler have also used an an gle- constraint action in their search for a massless 



phase in the 2-d 0(3) model 13j, Il4| . while Hasenbusch used the same action to ar- 



gue that M.P(N — 1) models are in the same universality class as O(N) models 15 



Refs. 13Ml5| presented numerical evidence that the action with the angle constraint 
falls in the same universality class as the standard action. Our study will confirm 
these results and will extend them by studying the cut-off effects of this topological 
action, as well as by investigating the topological susceptibility and the correlator 
of the topological charge density. Lattice actions with a similar constraint have also 



been used in |16l-l24|. however, not with an emphasis on the topological properties 
of some of these actions. The second topological lattice action that we consider re- 
ceives local contributions from the absolute value of the topological charge density. 
This action does not have the correct classical continuum limit either, but it obeys 
a lattice Schwarz inequality. Irrespective of this, as we will see, the correct quantum 
continuum limit is reached for both topological actions. 

It should be pointed out that lattice theories with topological actions are not 
regularizations of the topological field theories that arise in the context of string 



theory or conformal field theory [25] . While those theories realize new universality 
classes, the theories with topological lattice actions studied here fall into standard 
universality classes, despite the fact that they violate basic principles of classical 
physics. 

The O(N) model in (N — 1) dimensions has a non-trivial topological charge 
Q E TIn^IS"- 1 } = Z. We will investigate the 1-d 0(2) and the 2-d 0(3) model 
which have topological charges in Il^S 1 ] and n 2 [5 2 ], respectively. While the 1-d 
0(2) model will be studied analytically, the 2-d 0(3) model is investigated using 
Monte Carlo simulations. The 2-d 0(3) model can also be viewed as the N = 2 
member of the 2-d CP(N — 1) model family [26]. The CP(N — 1) manifold is the 
coset space SU(N)/U(N - 1) = S^^/S 1 . Since ^[S™- 1 ] = {0}, one obtains 
n 2 [CP(iV - 1)] = Tl^S 1 ] = Z, i.e. all 2-d CP(N - 1) models possess a non-trivial 
topological charge. Since they are asymptotically free, have an anomalously broken 
classical scale invariance and a dynamically generated mass gap, as well as instantons 
and ^-vacuum states, 2-d CP(N — 1) models share many features with 4-d Yang- 
Mills theories. This has motivated their detailed study beyond perturbation theory. 
The ^-vacuum physics of 4-d SU(N) gauge theories and of 2-d CP(N — 1) models 



has been reviewed in 27 . 



Topological aspects of lattice CP(N — 1) models have been investigated in 
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[si, B, E3, 28-35|. Depending on the lattice action and the lattice definition of the 
topological charge, the quantum continuum limit of the topological susceptibility 
Xt = (Q 2 ) /V, where V is the space-time volume, may be spoiled by short- distance 
lattice artifacts. These so-called dislocations have topological charge \Q\ = 1 and a 
minimal value of the lattice action. Semi-classical arguments, which are, however, 
not rigorous, suggest that Xt ma Y have a power-law divergence in the quantum con- 
tinuum limit |29|. This is expected to happen in the CP (2) model, when one uses 
the standard action in combination with the geometric definition of the topological 
charge [32)]. This problem does not arise for CP(N — 1) models with N > A. Even in 
the CP(2) case, dislocations can be eliminated and the proper continuum limit of Xt 
can be attained if one uses a modified lattice action [32] . Dislocations have also been 
identified in 4-d lattice Yang- Mills theory 36-38)]. Again, semi-classical arguments 
suggest that they may spoil the quantum continuum limit in SU{2) and 577(3) 
Yang-Mills theory, if the standard Wilson action is used in combination with the ge- 
ometric definition of the topological charge 0,0, BBS. As in the CP (2) model, 



dislocations can be eliminated in SU (2) and SU (3) Yang-Mills theory by using an 



improved lattice action |38J. The situation is more subtle in the 2-d 0(3) (or equiv- 
alently CP(1)) model. In this case, a semi-classical calculation in the continuum 
already yields a divergent topological susceptibility Xt even in a small space-time 
volume [30]. This divergence is not caused by lattice artifacts, but is an intrinsic 
feature of the theory in the continuum limit. Hence, one concludes that a meaningful 
quantum continuum limit of Xt does not exist in the 2-d 0(3) model. This is sup- 
ported by a calculation of Xt using a classically perfect lattice action in combination 



with a classically perfect topological charge |33J, which eliminates dislocations and 
indeed shows no power-like divergence of the topological susceptibility. However, 
Xt still diverges logarithmically, and thus a meaningful quantum continuum limit 
is not reached for this quantity. As we will see, the logarithmic divergence of Xt 
even arises for topological actions, although in that case one might have expected 
a power-law divergence due to dislocations. On the other hand, the correlator of 
the topological charge density will turn out to have a finite continuum limit. The 
concepts of classical and even quantum perfect definitions of the topological charge 



have also been investigated analytically in the 1-d 0(2) model 42 



In QCD, Ginsparg- Wilson lattice quarks |43| obey an Atiyah-Singer index theo 



rem even at finite lattice spacing 44]. Based on Ginsparg- Wilson lattice fermions, 
unambiguous definitions of the topological susceptibility which are free of short- 
distance singularities have been provided 45 f [46[. They have been used in a deriva- 
tion of the Witten-Veneziano formula l47H49j for the //-meson mass in a fully regu- 




larized non-perturbative framework 

This paper is organized as follows. Section 2 contains an analytic investiga- 
tion of the 1-d 0(2) model using two different topological lattice actions: one that 
suppresses topological charges and one that does not. In both cases, the correct 
quantum continuum limit is obtained. In Section 3 we study the 1-d 0(3) model in 
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a similar manner. Section 4 presents a Monte Carlo study of the 2-d 0(3) model. 
Again, we use two different topological lattice actions: one that does and one that 
does not obey a Schwarz inequality. As before, the correct quantum continuum limit 
is reached in both cases. The Monte Carlo data for the topological susceptibility 
are consistent with a logarithmic divergence, while the correlator of the topological 
charge density has a finite continuum limit. We summarize our results and draw 
conclusions in Section 5. 



In this section we consider the 1-d 0(2) model as an analytically solvable test case 
with two different topological lattice actions: one that explicitly suppresses topolog- 
ical charges and one that does not. Remarkably, irrespective of this, and although 
the lattice theories do not yield the correct classical limit, they do have the correct 
quantum continuum limit. 



In this subsection, we analytically solve the 1-d 0(2) model in the continuum. The 
results will then be compared with those of the corresponding lattice models. The 
1-d 0(2) model is equivalent to a quantum mechanical rotor. Let us consider a 
particle of mass M on a circle of radius R, and thus with the moment of inertia 
I = MR 2 . The Hamilton operator takes the form 



where <p is the angle describing the position of the particle, and 9 is analogous to the 
vacuum angle in QCD. At finite temperature T = 1//3, the corresponding Euclidean 
continuum action is given by 



2 The 1-d 0(2) Model 



2.1 The 1-d 0(2) Model in the Continuum 




(2.1) 




(2.2) 



where the topological charge takes the form 




P 



(2.3) 



The energy eigenfunctions of the Hamilton operator H{9) are given by 




(2.4) 
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where m G Z specifies the angular momentum, and the corresponding energy eigen- 
values are 

If 6 N :l 



The canonical partition function takes the form 

Z{6) = Tiexp(-f3H(9)) = ]T exp(-/9£ m (0)). (2.6) 

m<=Z 

The (not yet normalized) distribution of the topological charge Q is obtained as a 
Fourier transform of Z{6) 



p(Q) = i- jT ^ Z(0)exp(-i0Q) = ^ ^exp ("^V) , (2-7) 
and the topological susceptibility (evaluated at 9 = 0) reads 

x '- — - p T. Qa P(Q) ' ( ' 

In the zero temperature limit j3 — > oo this expression reduces to 

* = 3^7' (2 ' 9) 

The correlation length £ (again evaluated at = 0) is determined by the gap between 
the ground state and the first excited state 

^ l( o)- go( or 2J - < 2 ' 10 > 

such that at zero temperature 

** = S5- (2 ' n) 

Indeed, in the 1-d 0(2) model the topological susceptibility is a quantity with a 
meaningful finite quantum continuum limit, which scales like the inverse correlation 
length. As we will discuss later, this is not the case in the 2-d 0(3) model. 

It is interesting to minimize the action (at 6 = 0) in a given topological charge 
sector. In the 1-d 0(2) model the minimizing configurations take the form 

¥>(*) = ¥>(0) + (2.12) 

and they have the action 2n 2 IQ 2 / (3. Consequently, the action and the topological 
charge of all configurations obey the inequality 

> *w (2 . 13) 



It should be noted that, unlike instantons in 4-d non-Abelian gauge theories and in 
2-d CP(N — 1) models, in the 1-d 0(2) model the topologically non-trivial minimal 
action configurations are not concentrated at an instant in Euclidean time, but are 
homogeneously distributed over time. In this sense, they do not deserve to be called 
instantons. 



2.2 A Topological Lattice Action without Topological Charge 
Suppression 

Let us now consider a 1-d lattice 0(2) model with spin variables ip t G] — ir, ir], i.e. an 
XY model, with zero action for nearest neighbor spins with \(ip t+a — ip t ) mod 27r| < 5 
and infinite action for \(<pt+ a — ft) m od 2ix\ > 5. Here a is the lattice spacing. The 
geometric definition of the topological charge is given by 

i 

with ((ft+a — ft) m od 2n g] — n, ir]. The partition function takes the form 

Z(9) = Tr T(9) N , (2.15) 
with (5 = Na. The transfer matrix T(9) has the elements 



(ip t \T(9)\ip t+a ) = exp y-^iVt+a ~ ft) mod 2rr J , (2.16) 

for \((pt+ a ~ ft) mod 2tt| < 5. For \((pt+ a — ft) mod 2n\ > 5, on the other hand, the 
transfer matrix elements vanish. 

The transfer matrix can be diagonalized by changing to a basis of angular mo- 
mentum eigenstates 

(m\T(9)\m') = -— ^ J dip t J dip t+a (m\ipt)(ft\T(9)\ip t+a )(ipt + a\m') 



1 



dip t I dip t+a exp I -i-^—(ip t+a - ft) mod 2n 



(2tt) 2 J_ w J^_ s \ 2tt 

x exp {-irrnpt + im'ip t+a ) 

= $ mm > ^ J dipexp ^-i-^-ip + imip^j . (2.17) 

Hence, the eigenvalues of the transfer matrix are given by 

1 f S , (.( 9\ \ sin((m - 0/27r)<5) 



exp(-aE m (9)) = J dip exp (m - <pj 



(m-9/2ir)ir 

7i 6 V 2tt V 



(2.18) 



7 



In the continuum limit, a — > 0, we thus obtain 



Here we have identified the moment of inertia as 

3a 

I = -. (2.20) 

It should be noted that exp(— aE (0)) = 5/ir, such that the ground state energy- 
diverges in the continuum limit. This is no problem because only energy differ- 
ences are physically relevant (in this context). In order to reach finite results in the 
continuum limit a — > 0, we must put 5 = ^3a/I. Interestingly, in this limit, the 
topological lattice model reproduces the continuum 1-d 0(2) model. However, it 
should be noted that the lattice transfer matrix is not positive definite. In partic- 
ular, for \m — 8/2it\5 > tt /2 the transfer matrix eigenvalue exp(— aE m (8)) becomes 
negative. This is not necessarily problematical, as long as this behavior does not 
affect the continuum limit. The transfer matrix eigenvalues which obey 



9 

m 

2vr 



7T 7T I . 

< a = iV& (2 - 21 » 



are positive. This condition is automatically satisfied in the continuum limit a — > 0. 
We thus conclude that the topological model indeed provides an adequate regular- 
ization of the 1-d 0(2) model. 

It is interesting to investigate the cut-off effects of the topological lattice action. 
It is well-known that for the standard action the lattice artifacts set in at 0(a 2 ). In 
particular, the dimensionless ratio of energy gaps of the first two excited states is 
given by 

E s J0) - E°(0) ( a 2 a 3 \ 1 . 

' -4 1--3-3-3 + ... , 6= ? r\\ E^V . (2-22) 



Ef(0) - £ *(0) V 6 2 G, 8 7 ' £f(0) - E° Q (0) 

while for a classically perfect action the lattice artifacts are even exponentially sup- 
pressed [42] 



£f (0) - E%(0) \ n V vra ^ V 4a 7 J ' ^(q) - E§(0) ' 

(2.23) 

For the topological lattice action under consideration, the corresponding result takes 
the form 

«■(»)-«>(»)- 4 Cl + 55 + ...V (2.24) 



£,(())- £4(0) V 5f 



[42[ the corresponding expressions look different because they are expressed in terms of the 
continuum correlation length. 
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Because the topological lattice action does not obey the correct classical contin- 
uum limit, it suffers from strong lattice artifacts of 0(a). It should be noted that 
Symanzik's systematic effective theory for lattice artifacts 0, 0] is applicable in 
quantum field theory but not in quantum mechanics. Indeed, the terms of 0(a) in 
the topological action and of 0(a 3 ) in the standard action would be forbidden if 
Symanizik's theory would apply. The scaling behavior of the various lattice actions 
is illustrated in Figure 1. Indeed, one sees that the results obtained with the topo- 
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Figure 1: The scaling behavior of the ratio (E2 — Eq)/{E\ — Eq) of energy gaps 
for different lattice actions: standard action (solid curve), topological action with 
(dotted curve), and without topological charge suppression (dashed curve), as well as 
classically perfect action (dashed- dotted curve). The topological lattice actions suffer 
from cut-off effects of 0{a), while the standard action has only 0(a 2 ) artifacts. For 
the classically perfect action, the lattice artifacts are even exponentially suppressed. 
In the continuum (E 2 — E )/(E 1 — E ) = 4. 

logical lattice action converge much slower than the ones for the standard or the 
classically perfect action. Still, as pointed out before, the topological lattice action 
has the correct quantum continuum limit. 

Since the information about the topological charge distribution is encoded in the 
^-dependence of the energy spectrum E m (6), the topological susceptibility xt as well 
as other related topological quantities will also automatically come out correctly. In 
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order to show this explicitly, let us also consider the partition function of the lattice 
model 



Z{9) = Tr T(6) N -- 
such that 

p(Q) 



sin((m - 9/2n)5) 
(m-0/2n)7T 



N 



(2.25) 



— / d9 Z(d) exp(-i9Q) 



1 

1 

2^ 



de 



2sin(^/2vr) 



-i N 



exp( 



(2.26) 



Since p(Q) is the Fourier transform of an iV-th power, it is given by the iV-fold 
convolution p — pi -kpx ~k ■ ■ ■ *p\. The elementary distribution is given by 



Pl(Q) = 7T 
lix 



00 Jfl 2sin(gg/27r) 
at) exp( 

oo V 



5 



iOQ) = Q[Q + — - e Q - — 



2tt 



5 



2^ 



(2.27) 

In this case, Q e R is not restricted to integer values, because the elementary 
contributions to the total topological charge originate from local regions with open 
boundary conditions. We define the step function as Q(Q) = 1 for Q > 0, and 
as Q(Q) = otherwise. In the zero-temperature limit (3 — > oo, the topological 
susceptibility now takes the form 



Xt 



(2.28) 



which is indeed the correct result of eq. fl2.9p for the 1-d 0(2) model in the continuum. 



It is interesting to note that this lattice action violates the inequality ( I2.13P . This 
is obvious, because all allowed configurations have zero action. Still, for a given value 
of S and for a given inverse temperature /3 = Na, the allowed topological charges 
are restricted to 

\QM<^-N=^-, (2.29) 



such that 



2ir 2 IQ[<pf 



< 



2tt 
nI5 2 /3 



2ira 

37T/3 



— > oo. 



(2.30) 



p ~ 2a 2 2a 
Hence, unlike in the continuum theory, topologically non-trivial field configurations 
are not suppressed by the lattice action in the continuum limit a — > 0. Remarkably, 
nevertheless the lattice theory has the correct quantum continuum limit. 

As a further scaling test, we compare the cut-off effects of the topological sus- 
ceptibility in units of the mass gap, i.e. Xt£,, for the various lattice actions. For the 
standard action one obtains 



Xt £s 



2tt 2 



1 + 



a 



a 



(2.31) 



10 



while for the classically perfect action one finds 



XtCc 



2tt 2 



a 



7l z 



exp 



4a 



+ 



(2.32) 



For the topological action, on the other hand, we obtain 



4tt 2 3 



-i -l 



log 



sin 5 



2tt 2 



a 



(2.33) 



As before the lattice artifacts of the topological action are of 0(a), while they are 
of 0(a 2 ) for the standard action and exponentially suppressed for the classically 
perfect action. The results for the various actions are illustrated in Figure 2. 
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Figure 2: The scaling behavior of the product Xt£, °f the topological susceptibility and 
the correlation length for different lattice actions: standard action (solid curve), 
topological action with (dotted curve) and without topological charge suppression 
(dashed curve), as well as classically perfect action (dashed- dotted curve). The topo- 
logical lattice actions suffer from cut-off effects of 0(a), while the standard action 
has only 0(a 2 ) artifacts. For the classically perfect action, the lattice artifacts are 
exponentially suppressed. The continuum value amounts to XtC, — l/X^ 2 )- 
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2.3 A Topological Lattice Action with Topological Charge 
Suppression 



Let us now consider a topological lattice action which receives local contributions 
from the topological charge density, i. e. 



S [f\ = ^/MPt+a - Vt) mod 2n\. 



(2.34) 



A (dimensionless) coupling constant A > suppresses configurations with non-zero 
topological charge. In particular, by construction this action obeys the inequality 



S[<p) >2ttA|QM|, 



(2.35) 



which is inconsistent with the corresponding inequality (I2.13P of the continuum 
theory. 

In this case, the transfer matrix takes the form 



(<p t \T(0)\tp t+a ) = exp(-\\(ip t+a - ip t ) mod 27r|)exp ( -i—(ip t+a - tp t ) mod 2tt j . 

(2.36) 

It is again diagonalized in the basis of angular momentum eigenstates and the cor- 
responding eigenvalues are given by 



exp(-a£' m (6')) 



1 

2^ 
1 



d(pexp(— \\<p\) exp ( i ( m — — j <p 



x 



7rA 2 + (m-#/27r) 2 

A cos mil 

2 



|a — exp(— 7rA) 



m sin mir 

2tt / V 2 



For large A, i.e. in the continuum limit, one then obtains 



E m (9)-E (0) = -\og 
a 



1 + A^ m "27 



If 9 V 

— > —— m 

a\ 2 I 2tt/ 



In order to match the continuum result of eq. fl2.5p . we thus identify 

, a\- x 121 
I = — => A = 




.(2.37) 



(2.38) 



(2.39) 



Using the inequality (I2.35p . one then concludes that the action of topologically non- 
trivial field configurations diverges in the continuum limit, 



S[(p]>2n\\Q[(p]\ = 2icJ^\Q[(p]\ ->• oo. 



(2.40) 
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Remarkably, despite this fact, in the quantum continuum limit the ^-dependent 
energy spectrum still agrees with the one of the continuum theory. 



As before, we consider the lattice artifacts of the ratio of energy gaps, which now 
takes the form 

Again, the lattice artifacts, which are also illustrated in Figure 1, are of O(a). The 
artifacts of the topological lattice action with topological charge suppression are 
even a factor 5/2 larger than for the topological action without topological charge 
suppression. In particular, even for £ = 15a the deviations from the continuum limit 
are as large as 10 percent, while they are only about 0.5 percent with the standard 
action. 

Let us again consider the ^-dependent partition function 
Z{9) = TiT(9) N = J2^M-NaE m (9)) = 



mez 



Ell 



1 



TT A 2 + (m - 9/2ir) 2 



{A — exp(— 7rA) 



x 



A cos mn I — I m I sin I mn — 

2 V 2nJ V 2 



N 



(2.42) 



As before, the corresponding topological charge distribution p(Q) is obtained as the 
iV-fold convolution of an elementary distribution 



pi(Q) 



A — exp(— 7rA) 



9 9 9' 
A cos sin - 

2 TT 2 



= [1 - exp(-2vrA)] exp(-27rA|Q|). 
The topological susceptibility in the zero-temperature limit is then given by 



exp(— i9Q) 
(2.43) 



Xt = 



i 



af^dQ Pl (Q) 2n 2 \ 2 a 4% 2 f 
which again is the correct quantum continuum limit. 



(2.44) 



Let us again consider the lattice artifacts of the product XtC Up to exponentially 
small corrections, we obtain 



1 



2tt 2 A 2 



l0 § ( 1 + p 



2tt 2 



(2.45) 



As for the topological action of the previous subsection, the lattice artifacts are of 
0(a). This result is also illustrated in Figure 2. 
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We hence conclude that, despite the fact that the two topological lattice actions 
do not have the correct classical continuum limit, cannot be treated perturbatively, 
or violate the classical inequality ( I2.13p . they both have the correct quantum con- 
tinuum limit. In particular, this holds for the ^-dependent energy spectrum and all 
quantities derived from it, including the topological susceptibility. 



3 The 1-d 0(3) Model 

In this section we consider an angle-constraint topological lattice action for the 1-d 
0(3) model, which describes a quantum mechanical particle moving on the surface 
of a sphere S 2 . Again, despite the fact that the topological lattice action does 
not obey the correct classical continuum limit, it correctly reproduces the quantum 
continuum limit. 



3.1 A Particle Moving on S 2 

Let us now consider a particle of mass M moving on a sphere S 2 of radius R. The 
Hamiltonian then takes the form 

H-£, (3.1) 

where L is the angular momentum and I = MR 2 is the moment of inertia. The 
corresponding eigenfunctions are the spherical harmonics (9, ip\lm) = Yi m (6, (p) with 
/ G {0, 1,2,...} and m G {—I, — I + 1, . . . , I}, and the {21 + l)-fold degenerate energy 
eigenvalues are 

* - (3.2) 

In this case, one cannot construct a topological charge and the Euclidean action is 
simply given by 

S[e\ = I dt -d t e-d t e, (3.3) 



2 

with e(t) = (sin6'(t) cos ip (t), sin 6{t) siia ip(t), cos 9 (t)). 



3.2 The 1-d 0(3) Model with a Topological Lattice Action 

Let us now consider the 0(3) model with spins e t = (sm9 t cos<p t ,sm9 t sm(p t ,cos9 t ) 
attached to the sites t of a 1-d lattice with spacing a. The lattice action constrains the 
angle between neighboring spins e t and e t + a to a maximal value 5, i. e. e t -e t + a > cos S. 
The action vanishes as long as this constraint is satisfied and is infinite otherwise. 
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The corresponding transfer matrix is then given by (9 t ,(p t \T\9 t+a ,<p t+a ) = 1 for 
e* • e t+a > cos 5, and (9 t , tp t \T\B t+a , <p t+a ) = otherwise. We now put e t = (0,0, 1) 
and e t+a = (sin 6 cosy?, sin 9 sin <p, cos 6). Inserting complete sets of states \lm) and 
using the fact that the transfer matrix is 0(3)-invariant, i.e. 



(lm\T\l'm) = Sii>S m m' exp(-a£j), 

one obtains 

9(cos0-cos<5) = (0,0\T\e,<p) = ^{0,0\lm){lm\T\lm){lm\6,(p) 

l,m 

= J2 Y i™( >°) e M-aEi)Y lm (9,<py. 



(3.4) 



(3.5) 



Lin 



Inverting this relation, we find 
^ m (0,0)exp(-a£,) = 



1 p2n 

dcos9 / dp 6(cos# — cos 5)Yi m (9, tp) 
i Jo 



dx 



cos <5 



2Z + 1 



47T 



(3.6) 



where Pi(x) is a Legendre polynomial. Using Y/ m (0,0) = 5 m ,o^ (2/ + l)/47r and 
applying the Legendre differential equation (with P((x) = dPi{x)/dx) 



d_ 

dx 



[(i - x 2 )p;(x)} +i(i + i)Pi(x) = o, 



(3.7) 



for / 7^ one obtains 



exp(— aEi) = 2n dx Pi(x) 

J cos <5 



2tt 



1(1 + 1) 



2tt 



sin 2 5 P[ (cos 5). 



1(1 + 1) 

Similarly, for I — one finds exp(— aE ) = 2n(l — cos 5), which results in 



(3.8) 



Ei - E = — log 

a 



1 + cos 5 



Pl(cos5) 



1(1 + 1) 

Expanding in small values of 5, i.e. putting cos 5 ~ 1 — S 2 /2, one obtains 



(3.9) 



Ei — En 



--log 

a 



= — log 

a 



2 



1(1 + 1) 



p[(i)- u -[p!(i) + 2p;>(i)} 



+ 0(S 4 ) 



8a 



l(l + l) + 0(S 4 ). 



(3.10) 
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Here we have used 



^P/(l) + 2^(l)=(^) 2 



(3.11) 



Hence, identifying 




(3.12) 



we indeed reproduce the correct spectrum of eq. (13.21) in the quantum continuum 
limit. It turns out that in the 1-d O(N) model the moment of inertia is given by 



Let us again consider the lattice artifacts in the ratio of energy gaps, which now 
takes the form 



In this section we consider the 2-d 0(3) model using numerical simulations. We first 
summarize the results obtained before in the continuum and with standard as well 
as modified lattice actions. In analogy to Subsections 2.2 and 2.3 we then investigate 
two topological lattice actions, one that does and one that does not obey a Schwarz 
inequality between action and topological charge. Remarkably, in both cases we 
will again obtain the correct quantum continuum limit. While it is well-known that 
in this model the topological susceptibility is logarithmically divergent, we will see 
that the correlator of the topological charge density has a finite continuum limit. 
Hence there are topological quantities that do have a finite continuum limit in this 
model. 

4.1 The 2-d 0(3) Model in the Continuum and with the 
Standard Lattice Action 

In this subsection we summarize results obtained before either in the continuum or 
using the standard lattice action for the 2-d 0(3) model. In the continuum, the 2-d 
0(3) model has the Euclidean action 



/ = (N + l)a/S 2 . 




(3.13) 



As for the 1-d 0(2) model, the lattice artifacts are of 0(a). 



4 The 2-d 0(3) Model 




(4.1) 
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Here e(x) G S* 2 is a 3-component unit-vector field defined at each point x in a 
2- dimensional Euclidean space-time, and the topological charge is given by 

QW\ = -^ j d2x £ ^e • (fye x d v e) G U 2 [S 2 } = Z. (4.2) 

It is straightforward to show that the following integral, which is non-negative by 
construction, takes the form 

I = J d 2 x (<9 M e ± £ M „e x d v e) 2 = Ag 2 S[e\ ± 16nQ[e\ > 0. (4.3) 

This immediately implies the Schwarz inequality 

S[e\ > pQ[e\\. (4.4) 

Field configurations which saturate this inequality are (anti-)self-dual, i.e. 

d^e = ie^e x d v e, (4.5) 

and are known as (anti-)instantons. For these configurations the Lagrangian £(e) 
(at 6 = 0) is proportional to the absolute value of the topological charge density 
q(e), i.e. 

1 1 , „ | 47T 

= Y^ d ^' d ^ = 2g^\ £ ^' ( d ^ x 9 v e)\ = — \q{e)\. (4.6) 

In Section 4.3 we will introduce a topological lattice action such that C(e) is pro- 
portional to \q{e)\ for all configurations, not just for instantons or anti-instantons. 



Remarkably, at 9 = the 2-d 0(3) model can be solved exactly using the Bethe 



53H55| . Based on these results, using the Wiener-Hopf technique, the exact 



ansatz 

mass gap of the 2-d 0(3) model 

g 

m = -Ajgs, (4.7) 

has been derived [56[ . Here e is the base of the natural logarithm, and is the 
scale generated by dimensional transmutation in the modified minimal subtraction 
renormalization scheme. Even the finite-size effects of the mass gap m(L) for the 
2-d 0(3) model on a finite periodic spatial interval of size L have been calculated 
analytically A dimensionless physical quantity 

u = Lm(L) = L/£(L), (4.8) 

is then obtained as the ratio of the spatial size L and the finite-volume correlation 
length = l/m(L). Similarly, one defines the step scaling function (with scale 
factor 2) 

a(2,u ) = 2Lm(2L), (4.9) 



17 



which is thus also known analytically. Later, we will compare Monte Carlo data for 
cr(2,uo), obtained with a topological lattice action, with the analytic result. The 
step scaling function was first studied on the lattice by Liischer, Weisz, and Wolff 
[58J using the standard action 

Now x denotes the lattice sites, and p, is a vector of length a pointing in the p- 
direction. The step scaling function is affected by lattice artifacts which — for some 
time — seemed not to be described by Symanzik's effective theory. Interestingly, a 
recent careful study of the lattice artifacts has shown that large logarithms arise, 



and that Symanzik's theory does describe the lattice artifacts correctly [59|, |60 
Hence, one may conclude that the continuum limit of the finite-volume mass gap 
m(L) is finally well understood. Furthermore, excellent agreement between the 
Zamolodchikov bootstrap S-matrix and Monte Carlo data, obtained with both the 



standard and a classically perfect lattice action, has been reported [61 
Let us now discuss the topological susceptibility 

of the 2-d 0(3) model. Here V = f3L is the space-time volume. Based on naive 
dimensional analysis, one would expect that Xt£, 2 approaches a constant in the con- 
tinuum limit. While this expectation is met in 2-d CP(N — 1) models with N > 3, 
the CP(1) (or equivalently 0(3)) model behaves differently. The pathological behav- 
ior of Xt m the 2-d 0(3) model manifests itself already in the continuum formulation. 
As was pointed out by Liischer [29], the integration over the instanton size parameter 
p in a semi-classical calculation gives rise to a logarithmically divergent ultra-violet 
contribution to Xt which is proportional to J dp/ p. 

Berg and Liischer have investigated Xt using a geometric definition of the lattice 



topological charge 28]. In this definition, each plaquette of the square lattice is 
divided into two triangles, as illustrated in Figure 3. The spins e\, e 2 , and e*3 at the 
three corners of a lattice triangle ii 23 define the corners of a spherical triangle on 
S 2 . The oriented area A 123 of the spherical triangle is given by 

A X 2z = 2ip e [-27r,27r], X + iY = rexp(iip), 

X = 1 + ei ■ e 2 + e 2 • e 3 + e 3 • e u Y = e x ■ (e 2 x e 3 ). (4.12) 

The geometric topological charge is the sum of the oriented areas A xyz over all 
triangles t xyz , normalized by the area Air of S 2 , i.e. 
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Figure 3: Decomposition of the square lattice into triangles to be used to define the 
action and the topological charge. 

The decomposition of the square lattice into triangles illustrated in Figure 3 is in- 
variant under it/ 2 rotations and translation invariant by an even number of lattice 
spacings. Obviously, the topological charge inherits these symmetries. As illus- 
trated in Figure 4, an individual lattice plaquette 1234 can be divided into two 
triangles 123 and 243 or alternatively into triangles 124 and 143. In general, it is 
not guaranteed that A123 + A243 is equal to A124 + A143. However, if the relative 
angle between nearest-neighbor spins is smaller than 7r/2, one can indeed show that 
A 12 3 + A243 = A 12 4 + A U3 . In that case, the topological charge becomes indepen- 
dent of the particular decomposition into triangles, and thus becomes invariant even 
against translations by a single lattice spacing. Here we prefer to work on a triangu- 
lated quadratic (rather than a triangular) lattice because the definition of the step 
scaling function refers to a rectangular space-time volume. 

The geometric topological charge is undefined for a set of exceptional field config- 
urations, which form a set of measure zero in configuration space. These exceptional 
configurations contain a spherical triangle that covers exactly one half of S 2 and thus 
has an area ±2tt. The infinitesimal neighborhood of an exceptional configuration 
contains configurations whose topological charges differ by 1. If the relative angle 
between nearest-neighbor spins is smaller than ir/2, exceptional configurations can- 
not arise. When one uses the standard lattice action, exceptional configurations 
have a finite action. In the following subsection, we will consider a topological lat- 
tice action that constrains the relative angle between nearest-neighbor spins to a 
maximum value 5. When 5 < ir/2, exceptional configurations are excluded, and 
- just as in the continuum — different topological sectors are then separated by 
infinite-action barriers. This means that the angle-constraint topological action, but 
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1 



2 1 



2 



Figure 4: A given plaquette 1234 can be decomposed into two triangles in two alter- 
native ways: a) 123 and 243 or b) 124 an d 143- 



not the standard action, naturally leads to a unique segmentation of configuration 
space into distinct topological sectors. 

In addition to the logarithmic ultra-violet divergence that is present in the contin- 
uum theory, in the 2-d lattice 0(3) model Xt is affected by additional short-distance 
artifacts. These so-called dislocations are minimal action field configurations with 
a non-zero topological charge. As such, they depend on both the definition of the 
lattice action and the definition of the lattice topological charge. Using the standard 
lattice action in combination with the geometric definition of the lattice topological 
charge, the dislocations are exceptional configurations with an action Sd = c/g 2 
with c ~ 6.69 < 47r [29] . Semi-classical arguments (which, however, are not rigor- 
ous) suggest that the topological susceptibility should scale as Xt oc exp(— Sd)- Due 
to asymptotic freedom, the correlation length scales as £/a oc exp(27r / g 2 ) , where 2n 
is the universal 1-loop coefficient of the /3-function. Hence, one expects a power-law 
divergence of the dimensionless combination 

Xti 2 oc exp(-S d ) exp ( — J = exp [ g2 J ( ~ J ■ ( 4 - 14 ) 



Dislocations have been eliminated in [33| by using a classically perfect lattice action 
in combination with a classically perfect definition of the lattice topological charge. 
Perfect discretizations are based on the renormalization group and eliminate cut-off 
effects at the classical level. In particular, one then has Sd = 47r/g 2 , i.e. c = Air. 
Interestingly, the topological susceptibility was then found to still diverge logarith- 
mically. From all this one concludes that Xt does not have a finite continuum limit 
in the 2-d 0(3) model. 
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L/a 


5/n 


£(L)/a 


£{2L)/a 


10 


0.44858728 


9.4383(1) 


15.9515(15) 


16 


0.434009 


15.1013(2) 


25.4787(10) 


32 


0.415095 


30.2029(2) 


50.847(3) 


64 


0.398665 


60.406(1) 


101.578(4) 



Table 1: Constraint angle S that leads to uq = Lm(L) = 1.0595 for various lattice 
sizes together with the correlation lengths £(L) and £(2L). 



The situation is different in 2-d CP(N — 1) models with N > 3. While Xt 
is still divergent in the CP(2) model when one uses the standard lattice action in 
combination with the geometric definition of the topolo gica l charge, dislocations can 
be suppressed by choosing a modified lattice action 29|, |32|. For N > A, dislocations 
do not cause any problems, even when the standard action is used. Similarly, in 4-d 
SU(N) lattice Yang-Mills theories, when one uses the geometric topological charge, 
Xt suffers from dislocations for N = 2 and 3, which again can be suppressed by using 



a modified lattice action [37], [38 



4.2 A Topological Lattice Action without Topological Charge 
Suppression 

In analogy to the 1-d 0(2) and 0(3) models, we now investigate the 2-d 0(3) model 
with a topological lattice action that constrains the relative angle of nearest-neighbor 
spins on a square lattice to a maximum angle 5. All configurations that violate this 
constraint are forbidden and thus have infinite action, while all other configurations 
are allowed and have zero action. This action can be simulated with the very efficient 



Wolff cluster algorithm [62|, [63J]. Two neighboring spins are put in the same cluster 
if flipping one of them on a randomly chosen reflection plane would increase their 
relative angle beyond 5. For the efficiency of the algorithm it is essential that, using 
this method, only spins on the same side of the reflection plane end up in the same 
cluster. The correlation function of two spins (e x ■ e y ) can then be computed very 
accurately using an improved estimator. 

We have fine-tuned the maximal angle 5 such that the finite-volume mass gap 
m(L) satisfies uo = Lm(L) = 1.0595 for lattices with L/a = 10,16,32, and 64. 
The corresponding values of 5 are listed in Table 1. By measuring the mass gap 
m(2L) for these values, we have then determined the lattice value £(2,w ,a/L) = 
2Lm(2L) of the step scaling function, which is known to approach the continuum 



limit a(2,u = 1.0595) = 1.26121035(2) [59|, [60j. Figure 5 compares the cut-off 
effects of E(2, uq, a/ L) for the topological action with the standard and with two 
modified actions. At first glance, it seems that all four actions have lattice artifacts of 
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Figure 5: Cut-off dependence of the step scaling function S(2, tt , a/ L) atu = 1.0595 
for four different lattice actions: the standard as well as two different modified ac- 
tions (0(3)D(l/3) and 0(3)D(— 1/4) ) f^j /. and the topological lattice action without 
topological charge suppression. The lines are fits based on Symanzik's effective the- 
ory with the continuum value fixed to the exact result S(2, uq = 1.0595, a/L — > 0) = 
a(2,w = 1-0595) = 1.26121035(2). 



0(a). This would contradict Symanzik's effective theory. However, as investigated 



in detail |59|, l60j . the standard and the modified actions indeed have lattice arti- 



facts of 0(a 2 log 3 (a/L)), with large logarithmic corrections mimicking 0(a) effects. 
In particular, Symanzik's theory correctly describes the observed lattice artifacts. 
Since lattice perturbation theory is not applicable to topological lattice actions, one 
might think that one cannot use Symanzik's theory to predict the lattice artifacts. 
However, this is true only to some extent. Symanzik's effective theory is formulated 
in the continuum and should be applicable to any lattice theory that reaches the 
correct quantum continuum limit. As we have seen analytically for the 1-d 0(2) 
and 0(3) models, the topological lattice actions suffer from lattice artifacts of 0(a), 
which would contradict Symanzik's effective theory. However, since the underlying 
power-counting does not work in quantum mechanics, Symanzik's theory is not ap- 
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plicable in that case. In the 2-d case, however, Symanzik's theory applies even to 
topological lattice actions. This suggests to fit the lattice step scaling function to 

a 2 

S(2,m = 1.0595, a/L) = a(2,u = 1.0595) + — [B log 3 {L / a) + C log 2 {L / a) + ...], 

(4.15) 

which gives a good fit for B = —0.067(4) and C = 0.014(9). Interestingly, in the 
range considered here, the lattice artifacts of the topological action are smaller than 
those of the standard action, for which one obtains B = 0.041(6) and C = 0.26(2). 
In fact, for the standard action at L/a = 64 the sub- leading term proportional 
to log 2 (L/a) is still larger than the leading term proportional to log 3 (L/a), while 
this is not the case for the topological action. One can estimate that the lattice 
artifacts of the standard action will be smaller than the ones of the topological 
action only for correlation lengths larger than about 5 x 10 4 a. Furthermore, if 
one uses a(2,u = 1.0595) as a fit parameter, only the topological action data are 
consistent with the exact value within error bars, while the other actions give rise 
to small deviations. Given the fact that the topological action violates the classical 
continuum limit, and is thus tree-level impaired, it performs remarkably well. This 
may perhaps encourage the use of topological lattice actions also in other models, 
including Abelian and non-Abelian gauge theories. 

The very accurate approach to the exact continuum result for cr(2,M = 1.0595) 
strongly suggests that the topological lattice action indeed leads into the standard 



universality class of the 2-d 0(3) model. This confirms earlier results of [ill . ll3Ml5l 
and also justifies a posteriori the use of a topological lattice action in [1 



4.3 A Topological Lattice Action with Topological Charge 
Suppression 

In analogy to the action for the 1-d 0(2) model discussed in Section 2.3, we now 
introduce a topological lattice action for the 2-d 0(3) model which explicitly sup- 
presses topological charges and obeys a Schwarz inequality. Again, the action is 
given by the absolute value of the topological charge density, i. e. 

S[e\ = \J2\A xyz \. (4.16) 

txyz 

Here \A xyz \ is the area of the spherical triangle on S 2 defined by the spins e x , e y , 
and e z at the three corners of a lattice triangle t xyz , cf. eq. fl4.12l) . and A is a positive 
coupling constant. By construction, this action obeys the inequality 

S[e\ > \\Y,A*vz\ =47rA|Q[e]|. (4.17) 

txyz 

A comparison with the Schwarz inequality eq. (14.41) of the continuum theory may 
suggest to identify A = 1/g 2 , but, as we will see, this is not necessarily justified. By 
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construction, for this topological lattice action the Lagrangian C(e) is proportional 
to the absolute value of the topological charge density, i.e. C(e) = 47rA|g(e)|, for 
all configurations. In the continuum theory the corresponding eq. fl4.6p is satisfied 
only for instantons or anti-instantons. It is interesting to investigate the limit A — > 
oo. Then the allowed configurations only contain spherical triangles of zero area. 
Consequently, at A = oo all spins e x fall on a common great circle in S 2 , and thus 
seem to represent an 0(2) model. When a triangle t xyz contains an 0(2) vortex, 
the corresponding area is \A xyz \ = 2n. Hence, at A = oo, the 2-d 0(3) model from 
above reduces to a 2-d 0(2) model from which vortices have been eliminated. Such 
a model is expected to be in a massless phase. In the 2-d 0(3) model, the continuum 
limit is approached at A — > oo, not by putting A = oo. In particular, universality 
suggests that we should still recover the asymptotically free continuum limit of the 
2-d 0(3) model. 

Unfortunately, the action of eq. (I4.16P cannot be simulated with an efficient Wolff- 
type embedding cluster algorithm. While it is possible to define a cluster algorithm 
that is ergodic and obeys detailed balance, one is forced to put spins in one common 
cluster although they are on different sides of the reflection plane. This renders 
the algorithm inefficient. Hence we have used a Metropolis algorithm to simulate 
this action. A high-precision study of the mass gap, as described in the previous 
subsection, is then not feasible. A numerically better accessible quantity is the 
second moment correlation length, which does not require a fit of the correlation 
function G(x — y) — (e x ■ e y ) at large distances. Instead one considers the Fourier 
transform 

G(p) = G ( x ) exp(ipx) (4.18) 

X 

for a quadratic L x L lattice, which yields the susceptibility x — G(p = 0) as well 
as the corresponding quantity at the smallest non-zero momentum F = G(p = 
(27r/L,0)). The second moment correlation length is then defined as 

UL)=( A *~ F lT X 12 ■ (4.19) 



4Fsin 2 (7r/L) 

Using the corresponding mass rri2(L) = l/^L) one can define the step scaling func- 
tion for the second moment correlation length as £2(2, m,2(L)L, a/L) = 2Lm2(2L). 
In the continuum limit, ^{L)/a — > oo, this function is again universal and has 
been determined in The universality has been verified in the (2 + l)-d spin 

1/2 quantum Heisenberg model j66[, which dimensionally reduces to the 2-d 0(3) 



model in the low-temperature limit [67|. We have measured the step scaling func- 
tion £2(2, rri2{L)L, a/L) on quadratic lattices with L/a = 50 and 100. In addition, 
we have measured this quantity for the topological lattice action of the previous 
subsection, which does not explicitly suppress topological charges. As illustrated in 
Figure 6, in both cases one finds very good agreement with the step scaling function 
obtained with the standard action. This again confirms that topological lattice ac- 
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tions lead to the correct quantum continuum limit, despite the fact that the classical 
continuum limit is not correctly represented. 




Figure 6: The ratio of second moment correlation lengths ^(^L) / ^{L) obtained 
with the topological lattice action with (triangles) and without topological charge 
suppression (circles). Within error bars the data fall onto the universal curve that 
was extracted from simulations with the standard action WaJ . 



4.4 Topological Susceptibility from Topological Lattice Ac- 
tions 

As discussed in the introduction, in the 2-d 0(3) model the topological susceptibility 
Xt does not have a finite continuum limit. When the standard lattice action is used, 
based on (non-rigorous) semi-classical arguments one would expect a power-law di- 
vergence of Xt due to dislocations which are short-range lattice artifacts carrying 
non-zero topological charge j29[. Even when a classically perfect action is used in 
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combination with a classically perfect topological charge, Xt still diverges logarith- 



mically 33]. The logarithmic divergence is not a lattice artifact, but is an inherent 
feature of the 2-d 0(3) model even in the continuum. In this subsection we investi- 
gate the topological susceptibility using the two topological lattice actions with and 
without explicit topological charge suppression. 

First, we consider the topological lattice action of Subsection 4.2 which does not 
explicitly suppress topological charges and which violates the Schwarz inequality. In 
order to hold the physical volume fixed while approaching the continuum limit, we 
demand L = A&iL). Since it is easier to compute numerically, we base this criterion 
on the second moment correlation length ^{L), and not on the inverse mass gap 
£(L), which is consistently just about 1.7 percent larger than ^(L). In the infinite 
volume limit the ratio £/£2 = 1.0007(1) has been determined very accurately in 
j64| . As illustrated in Figure 7, the discrepancy between this result and the ratio 
£(Z/)/£ 2 C£) = 1.017(3) observed at L = 4£ 2 (L) can be attributed to finite volume 
effects. I 



1.02 



1.015 



1.01 



1.005 



0.995 




200 



Figure 7: The ratio ^(L)/^ 2 (L) as a function of L/a (at fixed 5/tt = 0.4849,) is 
consistent with the infinite volume result £/£ 2 = 1.0007(1) obtained in }64] - The 
infinite volume result is represented by the horizontal line. 



2 We thank the referee for bringing this issue to our attention. 
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Figure 8 shows the finite volume correlation length ^ 2 {L)/a as a function of the 
constraint angle 5, together with a fit of the form 

£ 2 (L) t 9 /2tt\ lb , 
' -Ag 2 exp( — ), - = - + c , (4.20) 



a \9 J 9 o 

which yields A = 0.24(2), b = 0.263(2), and c = -0.54(2), with xVd.o.f. « 0.7. 
Only the large lattices with L/a > 200 are included in the fit. The exponential 
increase of the correlation length is a consequence of asymptotic freedom. The fit 
assumes asymptotic scaling and contains the universal 1- and 2-loop coefficients 
of the /3-function. The relation between S and g is inspired by the corresponding 
relation eq. fl3.12p in the 1-d 0(3) model. It should, however, be pointed out that 
it is just a phenomenological ansatz which cannot be derived analytically, because 
perturbation theory does not apply to topological lattice actions. For the standard 
action, asymptotic scaling is known to set in only at very large correlation lengths 



651 ] . Hence, also the value of b determined here may not yet correspond to the 
asymptotic value in the continuum limit. Again, since perturbation theory is not 
applicable to topological lattice actions, one cannot determine the asymptotic value 
of b analytically. 

In order to investigate the scaling of the topological susceptibility, we consider 
the dimensionless combination 

16x t (L)6(L) 2 = = ( g 2)(L) , (4.21) 

as we approach the continuum limit ^2(L)/a — > oo. In Table 2 we list the re- 
sults of Monte Carlo simulations on various lattice sizes ranging from L/a = 40 to 
L/a = 800. Since this topological action vanishes for all allowed configurations, the 
dislocation action is Sd = 0. Hence, based on a naive semi-classical argument, one 
might expect \t^ 2 oc exp(— Sd)(£/a,) 2 = (£/a) 2 . Such a power-law divergence is not 
reflected by the numerical data depicted in Figure 9. Instead, the data are well 
fitted by 



(Q 2 >(L) = Alog^-^j, (4.22) 

with A = 0.91(1), L /a = 39(2), L 1 /a = 24(1), and x 2 /d-o.f. « 0.5. A power-law fit 
(Q 2 )(L) = BL U + C yields a small power v = 0.21(2), and x7d-o.f. w 1.9. 

A logarithmic divergence of (Q 2 )(L) was already encountered in the continuum 



theory [3CJ and was also observed on the lattice using a classically perfect action 
331 ] . It is interesting that the same behavior arises for the topological action which 
should be most vulnerable by dislocations. While the naive semi-classical argument 
is not rigorous, it is still remarkable that even the presence of zero-action dislocations 
does not spoil the logarithmic divergence of the continuum theory. This suggests 
that dislocations may not cause power-law divergences in other cases, including 2-d 
CP(N — 1) models and 4-d non-Abelian gauge theories, either. Indeed, no such 
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Figure 8: Second moment correlation length ^{L) /a at fixed physical size L = 4£ 2 (L) 
as a function of the constraint angle 5. The line is a fit of the form of eg. \J^.20^ , 
&(L)/a = Ag 2 exp( 2 ^/9 2 ) with 1/g 2 = b/5 2 + c. 



divergence has been detected, for example, in numerical data for Xt in 4-d 577(2) 



Yang-Mills theory [68 



We have investigated (Q 2 )(L) at fixed physical size L = 4£ 2 (£), also for the 
topological action of Subsection 4.3, which explicitly suppresses topological charges. 
The corresponding numerical results are listed in Table 3. Since in this case no 
efficient cluster algorithm is available, the calculation is limited to lattices up to 
L/a = 160. As illustrated in Figure 10, the data for (Q 2 )(L) are again consistent 
with the logarithmic divergence of eq. fl4.22p . Here the best fit yields A = 0.63(2), 
L /a = 20(3), Lx/a = 10(1), with x*/d.o.L w 0.4. 

In this case, we fit the second moment correlation length to the form 

M£> =Vexp (^\ J_ =6A , (4.23) 
a \g 2 J g 2 

which yields A = 0.19(1), b = 0.045(1). The results are illustrated in Figure 11. It 
should be pointed out that the relation between A and g is again just an ansatz, 
which cannot be derived analytically, because perturbation theory is not applicable 
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L/a 


5/tx 




i{L)/a 


(Q 2 )(L) 


40 


0.50405 


10.016(1) 


10.20(2) 


1.083(5) 


60 


0.48490 


14.996(8) 


15.29(4) 


1.299(6) 


80 


0.47260 


19.981(4) 


20.34(6) 


1.455(7) 


100 


0.46370 


24.99(2) 


25.46(8) 


1.598(9) 


120 


0.45680 


30.061(4) 


30.57(8) 


1.72(2) 


160 


0.44680 


39.98(2) 


40.67(9) 


1.923(9) 


200 


0.43950 


50.02(1) 


50.75(9) 


2.09(2) 


240 


0.43385 


59.98(2) 


60.97(9) 


2.23(1) 


320 


0.42545 


79.98(3) 


81.18(9) 


2.45(2) 


400 


0.41930 


100.07(3) 


101.5(2) 


2.64(2) 


480 


0.41455 


119.92(4) 


121.8(3) 


2.78(2) 


640 


0.40740 


159.87(7) 


162.5(3) 


3.046(8) 


800 


0.40210 


200.27(6) 


203.7(3) 


3.22(2) 



Table 2: Monte Carlo data for the topological charge squared (Q 2 )(L) at fixed physical 
size L = 4^2 (L) approaching the continuum limit ^2{L)/a — > oo. Here ^(L) is the 
second moment correlation length, defined in eq. \J^.19^ , while C,(L) is the inverse 
mass gap, and 5 is the constraint angle in the topological action of Subsection 4-2. 



L/a 


Atx\ 




(Q 2 )(L) 


40 


11.9215 


10.00(2) 


1.102(2) 


60 


13.781 


14.98(2) 


1.282(2) 


80 


15.112 


20.00(4) 


1.419(3) 


120 


16.988 


30.04(5) 


1.631(3) 


160 


18.325 


40.00(6) 


1.784(4) 



Table 3: Monte Carlo data for the topological charge squared (Q 2 )(L) at fixed physical 
size L = 4^2 (L) approaching the continuum limit ^2{L)/a — > oo. Here ^{L) is the 
second moment correlation length and A is the coupling constant in the topological 
action of Subsection 4-3. 
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L/a 



Figure 9: The topological charge squared (Q 2 )(L) at fixed physical size L = 4^2{L), 
for the topological lattice action of Subsection 4-2. As one approaches the continuum 
limit ^2(L)/a — > oo, (Q 2 )(L) diverges logarithmically. The line is a fit of the form 
ofeq.UM, (Q 2 )(L) = A\og((L + L )/L 1 ). 

to topological lattice actions. Furthermore, one cannot be sure that the above 
value of b persists in the continuum limit. The lattice Schwarz inequality (12.351) . 
S[e\ > 4ir\\Q[e\\, then translates into the inequality 

s ® > ( 4 - 24 ) 

which may be compared with the Schwarz inequality (14. 4p of the continuum theory 
S[e\ > (47i/g 2 )\Q[e\\. Since b ^ 1, the lattice Schwarz inequality deviates from 
the one of the classical continuum theory. However, the lattice theory still has the 
correct quantum continuum limit. 

4.5 Correlation Function of the Topological Charge Density 

Although the topological susceptibility Xt is logarithmically divergent in the 2-d 
0(3) model, this does not imply that the whole concept of topology is meaningless 



30 




200 



Figure 10: The topological charge squared (Q 2 )(L) at fixed physical size L = 4£ 2 {L), 
for the topological lattice action of Subsection 4-3. The line is a fit of the form of 
eq.UM, (Q 2 )(L)=Alog((L + L )/L 1 ). 



in the continuum limit. In particular, the correlation function (q{0)q(x)) of the 
topological charge density 



q(x) 



1 



57T 



[<9 M e(x) x d v e(x)] 



(4.25) 



69. 70 



Analytic 
The 



whose integral over x is Xt, has a finite continuum limit for \x\ > 0. 
results for this quantity have been derived by Balog and Niedermaier 
correlator is negative, except at x = 0. Up to logarithmic corrections, at short 
distances it has a power-law divergence proportional to \x\~~ A . At x = there is 
a positive divergent contact term. When the correlator is integrated over x, the 
numerical evidence obtained in the previous subsection suggests that the power- 
law divergence cancels against the contact term, but a logarithmic divergence of 



Xt persists. As shown analytically in [45|, |46j, using Ginsparg- Wilson quarks, also 
in QCD the corresponding short-distance power-law divergences cancel against the 
contact terms. A corresponding study in the large N limit of CP(N — 1) models 
has been presented in [34J . 



31 



" 30 




Figure 11: Second moment correlation length ^ 2 {L)/a at fixed physical size L 
4^2 (L) as a function of AttX. The line is a fit of the form of eg. lji4-23\ ), £ 2 (L)/a 
Ag 2 exp(27r / g 2 ) with 1/g 2 = bX. 



Let us consider the point-to-time-slice correlator 



G(x 2 ) 



dxi(q(0)q(x)) , x = (x 1 ,x 2 ). 



(4.26) 



The corresponding quantity on the lattice receives contributions A xyz /4ir from all 
triangles t Tyz in a row of plaquettes at fixed time x 2 . Using the meron-cluster 



algorithm [17j, we have constructed an improved estimator for (q(O)q(x)), which 
receives cluster-intrinsic contributions only. We have measured the correlator using 
the topological angle-constraint action with 5 /it = 0.4568 at L/a = 200, which 
yields £/a = 30.9(1), as well as with 5/tt = 0.4849 at L/a = 100, which yields 
£/a = 15.54(5). Assuming 0(a 2 ) cut-off effects, we have extrapolated these data 
to the continuum limit. In Figure 12, the extrapolated data are compared with the 
analytic results of 69, 70]U Finite size effects are expected to be small because 
we have worked at L f» 6£. At large distances, the extrapolated Monte Carlo data 
agree very well with the analytic prediction, while at short distances, x 2 j^ < 0.35, 
there are systematic deviations . We attribute these deviations to corrections to the 



3 We thank J. Balog for providing the numerical evaluation of the corresponding analytic results. 
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assumed 0(a 2 ) behavior, similar to the ones discussed before for the mass gap. In 
order to fully understand the cut-off effects, an analytic analysis along the lines of 



59l . |60| , combined with numerical data closer to the continuum limit would be most 
welcome. Keeping this in mind, we still conclude that our current data confirm 
again that the topological action leads into the standard 0(3) universality class. 
For the standard action a similar agreement had already been observed in 70]. In 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 12: The point-to -time- slice correlator G{x2) = J dxi(q(0)q(x)} of the topo- 
logical charge density (in units ofm 3 = l/£ 3 j as a function of Xijt,- The data points 
are extrapolations of Monte Carlo data obtained at£/a = 30.9(1) and^/a = 15.54(5) 
to the continuum limit, assuming 0(a 2 ) cut-off effects. At large distances, they are 
in good agreement with the analytic results of W f\ . \7U i represented by the solid curve. 
We attribute the systematic deviations at short distances to not fully understood 
cut-off effects. 



particular, this shows that some topological quantities make perfect sense in the 
continuum limit of the 2-d 0(3) model, despite the fact that Xt is logarithmically 
divergent. The divergence is expected to also affect the energy density of 9- vacua. 
Still, other physical quantities like the ^-dependent mass gap should have a well- 
defined continuum limit, which is accessible using the meron-cluster algorithm. 
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5 Conclusions 



We have investigated topological lattice actions for the 1-d 0(2) and 0(3) as well as 
for the 2-d 0(3) model. These actions are invariant against small deformations of 
the fields. Despite the fact that topological lattice actions do not have the correct 
classical continuum limit, as we have seen, they still yield the correct quantum 
continuum limit, irrespective of whether or not they explicitly suppress topological 
charges. In particular, it does not matter whether a topological action respects or 
violates a Schwarz inequality. Since, in contrast to other lattice actions, topological 
actions are invariant against small local deformations of the fields, one may have 
expected that they fall into a different universality class. However, the allowed 
local deformations are field-dependent and thus do not constitute a proper gauge 
symmetry of topological lattice models. In fact, even the standard action of a lattice 
0(N) model has a field- dependent local 0(N — 1) symmetry, since every spin can be 
rotated around the direction defined by the average of its nearest neighbors without 
changing the action value. Since such field-dependent local symmetries do not have 
the status of proper gauge symmetries, they have no impact on the corresponding 
universality class. 

Since topological lattice actions do not suppress small fluctuations of the fields, 
perturbation theory is not applicable. We have seen that in one dimension topo- 
logical lattice actions suffer from strong lattice artifacts of 0(a). This seems to 
contradict Symanzik's effective theory, which, however, does not apply in quantum 
mechanics. Despite the fact that lattice perturbation theory cannot be applied to 
topological lattice actions, Symanzik's effective theory, which is formulated in the 
continuum, still describes the lattice artifacts of the 2-d 0(3) model with a topo- 
logical lattice action. Interestingly, in contrast to the 1-d case, in the 2-d 0(3) 
model cut-off effects were observed to be less severe for a topological action than 
for the standard action, at least at practically accessible correlation lengths. Our 
results may encourage the use of unconventional regularizations, which might be 
advantageous from a computational point of view. 

Although the topological angle-constraint action does not explicitly suppress dis- 
locations, the corresponding topological susceptibility \t seems not to suffer from 
power-law divergences. Instead, the Monte Carlo data for Xt are consistent with a 
logarithmic divergence, which already occurs in the continuum theory. The numer- 
ical results for the point-to-time-slice correlator G(x 2 ) = L dx] (q(O)q(x)) of the 



topological charge density are consistent with the analytic predictions of |69|, [TO 
This underscores that, despite the fact that Xt is logarithmically divergent, there 
exist topological physical quantities that have a well-defined continuum limit in the 
2-d 0(3) model. Since the corresponding complex action problem can be solved 
using the meron-cluster algorithm, a numerical investigation of ^-vacuum effects in 
the 2-d 0(3) model is both feasible and physically meaningful. 
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It is straightforward to construct topological lattice actions for Abelian and non- 
Abelian gauge theories. One may simply constrain the trace of a Wilson plaquette 
variable by some minimal value. Configurations that satisfy this constraint on all 
plaquettes can then be assigned a zero action value. It is conceivable that one can 
take algorithmic advantages from actions of this kind. It is an interesting subject for 
future studies to decide whether topological lattice actions may be useful in lattice 
Yang-Mills theory or in lattice QCD. 

Our study underscores the robustness of universality, which, in particular, does 
not rely on classical concepts. Even when one uses actions that do not have the 
correct classical continuum limit, cannot be treated with perturbation theory, or do 
not obey a Schwarz inequality, the emerging quantum theory still has the correct 
continuum limit. This means that the standard approach of starting from a classical 
system and then quantizing it afterwards is not the only way to define a quantum 
theory. Even based on concepts that make no sense classically, one may still be 
able to construct a sensible quantum theory. Classical physics will then emerge 
dynamically from the underlying quantum system. 
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